Scalar mesons in rf — > 3tt°, n°7i + n decays 
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The decays rf —¥ and 7r°7r + 7r~ were consider in isobar model approach. Branching ratio 
and Dalitz-plot shape of the decay rf — > 3tt° can be explain by aO- and u-meson contributions. 
The a meson term is needed for the Dalitz-plot distribution reconstruction. The branching ratio of 
the decay rf — > tt°it + tv~ can be calculated in the same isobar framework. The prediction for the 
Dalitz-plot of this decay is shown which strongly dependents from model parameters. 
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I. INTRODUCTION 

Light scalar mesons are the subject of extensive theoretical and experimental investigation in recent years (see. [l[ 
and references therein). The structure of these mesons in terms of constituent quarks and gluons in rather vague. 
There are several models for light scalar mesons: usual quark- ant iquark states, tetraquark, meson-meson molecules, 
etc. Glueball states are also expected in w 1 GeV mass region scalar sector. 

The problem of light scalar mesons studying is that some of these states (a- and K-mesons) have large widths 
comparable with their mass, while the other (do (980) and /o(980)) lie near inelastic thresholds. For this reason one 
cannot expect usual Breit-Wigner form of the signal. There is also one class of reactions (for example rf —¥ rjim 
decay), that require the contributions of scalar mesons. Usually such reactions are considered in the framework of 
Chiral Perturbation theory 0-Q)- although for large energies the use of this formalism could be unacceptable. 

We think, that there is the other, more physically motivated approach to such decays. In the framework of so called 
isobar model the matrix element is saturated by contributions of virtual mesons with suitable masses and quantum 
numbers. In a series of works ([Ml) it was shown, that in order to describe the decay rf — > r/TTir it is sufficient 
to take into account the contributions of a- and ao-mesons. The branching fraction of this decay is reproduced by 
ao-meson contribution, and in order to explain the form of the Dalitz-plot, the contribution of a-meson should also 
be taken into consideration. 

In our article we use isobar model to study rf — >• 37r°, tt°tt + tt~ decays. It is well known that in these reactions G- 
parity is violated. According to Sutherland theorem [l3j the electromagnetic interaction in this case can be neglected, 
so only the difference of u- and d-quark masses can give the isospin violation. For this reason the study of rf — > 3tt 
decays could give valuable information on the ratios of current quark masses. In the framework of isobar model isospin 
violation is naturally parametrized in the terms of mixing of states with different isospin values and small mixing 
parameters should be proportional to m<j — m u mass difference. In our article we determine these parameters from 
fit of experimental branching fractions and Dalitz-plot distributions of rf — > 3ir°, ir°ir + Tr~ decays. 

The rest of the paper is organized as follows. In the next section we describe in detail the mechanism of isospin 
violation in the framework of isobar model and present diagrams that give contribution to rf — » 3ir°, 7r°7r + 7r~ decays. 
In sec. Ill and IV we use obtained matrix elements to analyze neutral and charged modes of rf — > 3ir decays. Short 
discussion is presented in Conclusion. In the Appendix we give the explicit form of used in our paper a-, do-, and 
p-meson exchange amplitudes. 

II. MODELS OF ISOSPIN VIOLATION 

Let us consider isospin violation decays rf — > 3n in the framework of isobar model, when the amplitude of the 
process is saturated by the contributions of virtual mesons with suitable masses and quantum numbers. 

This approach was used earlier to describe the isospin-conserving decay rf — > 777r°7r° In this works it 

was shown, that in order to reproduce the form of experimental Dalitz-plot it is sufficient to take into account the 
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FIG. 1: Doubly isospin violated diagrams of r/' 
tabular value 



3-7T decay. The particles shown with double line have isospin different from 



contributions of two scalar resonances: c-meson in 7r°7r°-channel and ao-meson in 7r°7y-channel. The matrix element 
of this decay can be written in the form 

— Aa (^7T7r) ~T" Aa (s?]Tri ) ~t~ Aa (^^2) ) 

where sy are invariant masses of corresponding pairs squared and A a (s) , A a (s) are the rescattering amplitudes in 7T7r- 
and ^-channels (explicit expressions for these amplitudes can be found in the Appendix) . fn our analysis of rf — > Stt 
decays we also restrict ourselves to contributions of these two resonances. 

When describing the isospin violation in considered here decays it is necessary to recall, that physically observed 
mesons have no definite value of isospin /. For example, ir°- and 77-mesons are actually superposition of / = and 
1=1 states 



phys 



K°>7=1 + e l 7r °)/=o' l r ?)phys = l»7>l=0 + 6 fo>J=l ■ 



The wave functions of other involved in the reaction particles can also be written in a similar form. The admixtures 
with "wrong" isospin values (| 7 r°) / _ and \rj) I=1 in the presented above examples) should be suppressed by the small 
factor 
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141. In what follows we 



that is usually interpreted as the mixing angle between states with different isospin values 
will consider only linear in this parameter terms. 

This restriction reduces significantly the number of diagrams that can give contribution to considered here processes. 
For, example, one can neglect the diagram shown in fig. [T^i, since in the decay | cr) z=0 — > | 7r °) J _ 1 | 7r °)/-o is° s P m is 
violated. As a result the corresponding diagram is quadratic in e and it can be neglected. The same is valid also for 
the diagram shown in fig. [TJd, where additional suppression factor stems from generalized Bose-Einstein symmetry, 
that in the isospin-conservation limit forbids the situation when two 7r-mesons are in the states with zero orbital 
momentum and total isospin I^n = 1- For this reason in the vertex |ao) J=1 — > I 71 " )/-! isospin is violated. 

Using these arguments it is easy to see, that only diagrams shown in fig. [2] can give linear in isospin violation 
parameter contributions to the considered in our article processes. The explicit form of the matrix element depends 
on what final state (37r° or 7r 7r + 7r~) we are talking about. In the case of neutral final state 



(1) 



the amplitude should be symmetric with respect to momenta of all final particles. In is clear, that it can be written 
in the form 



A [37T°] = e a A (si) + e a A a (si) + (si -> s 2 ) + ( Sl -> s 3 ) 



(2) 
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FIG. 2: Diagrams of 7/ — > 3tt decays, that violates isospin only once. The particles shown with double line have isospin 
different from tabular value 



where Sy = (pi + Pj) , amplitudes A a ,a{s) are given in the Appendix, and small parameters e a ,a will be determined 
from the fit of experimental data. In the case of charge final state 

7/ -> ir+ { Pi )tt- {p 2 )ir°{p 3 ) (3) 

no additional symmetrization is required. It is evident, that the particle with 1 = (e.g. c-meson on diagrams ®L,b 
or the upper of 7r-mesons on diagram [2}:) is neutral. So the amplitude of the decay ([3]) has the form 

A[tt°tt + tt-] = e a Av(s 12 ) + e a {Aa{si3)+Aa(s 2 3)}- (4) 

In some works (see [l5j]) p-meson contribution to these processes is also considered. It is clear, that the amplitude 
of = eta' — > 37T° decay cannot contain contribution from p-meson exchange, but in the case of charged final state this 
is not the fact. So, we should add a term e p A p (S12, S13) to the amplitude ((4]). The explicit expression for A p is given 
in the Appendix. 



III. 7/ ->■ 3tt° 

Let us first consider the neutral decay mode rj' — > 3ir°. We will approximate the matrix element of this decay by 
the following expression 

A(r]'^3ir ) = e a {Aa (su) + Aa (su) + Aa (S23)} + e a {A a {su) + A a {sis) + A a (S23)} ■ (5) 

Here Sjj = (pi + pj) 2 , A a (s), A a {s) are amplitudes of a- and arj-meson exchange (see Appendix for details), and e a , a 
are small isospin violation parameters that will be determined from experimental information about the branching 
fraction and dalitz-plot shape of this decay. 

The branching fraction of the decay r/ —> 37r° is equal to [l6| 

BR exp (?/ 3tt°) = (1.61 ± 0.23) x 1CT 3 . (6) 

The dependence of this branching fractions on the model parameters can be written in the form 

BRth = e 2 a B aa + e 2 a B aa + 2e cr e a B aa , (7) 

where B fftr , B aa and B aa are branching fractions of the considered here decays with only er-meson, ao-meson and a— cto 
interference contributions taken into account. Using presented in the Appendix matrix elements we get numerical 
values of these coefficients, presented in the second column of table HI 

The form of the dalitz-plot is usually written as an expansion in a variable Z defined as 

• ■ M n x 1 



(M v ' - 3to„.) 2 i= , 
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TABLE I: Contributions of different channels to the branching fraction and slope coefficient for rf 
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FIG. 3: Allowed parameters (e CT ,e a ) region for 77' — > 3ty° . The widths of the ellipse is equal to 2a (the error in branching 
fraction of this decay). The shaded arcs show the regions with given confidence level. 



Where Ti,2,3 are kinetic energies of 7r-mesons in initial meson rest frame. The parameter /3 in the expansion 

\M\ 2 = N(l + 2/3Z + ...). 
is called the slope coefficient of the dalitz-plot. Experimental value of this coefficient is [l|| 

/3exp (V -> 3tt°) = -0.59 ±0.18. 
The dependence of the slope coefficient on the model parameters can be written in the form 



2e a e a N a 



(8) 



(9) 



(10) 



where N aa , N aa , N aa and /3 CTCr , fi aa , (3 aa are the parameters of the expansions ([5} of squared matrix elements with only 
cr-meson, ao-meson and a — ao-intcrfcrcncc contributions taken into account. Numerical values of these parameters 
are presented in the third and fourth columns of table Q] It is clearly seen, that the slope coefficient of the dalitz-plot 
with only ao-contribution taken into account (i.e. parameter j3 aa ) is small. In other words, this resonance gives 
almost flat distribution over the dalitz plot, so it is necessary to consider also cr-meson to explain the experimental 
dependence of the squared matrix element on the invariant masses. Similar situation is observed also for 77' — > r]ir ir° 
decay, where for explanation of experimental dalitz plot cr-meson play crucial role [I2J]. It can also be noted, that the 
slope coefficient depends only on the ratio e a /e a . In fig. [4] we show the dependence of the slope coefficient on this 
ratio and demonstrate existing experimental constraints on /3. 

In order to have branching fraction if the decay 77' — > 3ir° consistent with the experimental value, parameters e a , a 
should lie in the region bounded by two ellipses shown in fig. [3] The width of this region is determined by the 
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FIG. 4: Dalitz plot slope coefficient of rj' 3ir versus the model parameter ratio 

experimental error of the decay ©■ It can be clearly seen, that by an order of magnitude this region coincides with 
the theoretical estimates presented in 14%. In our analysis we do not restrict ourselves to linear parametrization (|SJ) 



and fitted experimental dalitz-plot by the resonance parametrization ([5]) with fitting parameters e CTia . It turns out the 
at a = 0.1 confidence level 1 there are two symmetric regions on the ellipse (fig. [3]) where required confidence level can 
be achieved. 

These regions are rather large. On the one hand this is because the matrix element of the decay is not very sensitive 
to the ratio e a /e a . On the other hand, the large width of these regions reflects significant experimental errors in slope 
coefficient of the matrix element, caused by experimental difficulties of the considered decay. For example, in order 
to remove intensive background signals one has to remove some regions from the dalitz-plot [16[ . Moreover, originally 
the slope coefficient was extracted from experimental data assuming the the squared matrix element is linear in the 
variable Z, and our subsequent analysis shows the presence of higher non-linear terms. 

In fig. [3] the allowed regions of model parameter values are shown by shaded ellipse parts. As it was mentioned 
above, one can change these parameters without noticeable decrease in experimental data description. In order to 
demonstrate this fact and to check Z-distribution sensitivity to parameter values we considered two cases: equal 
values of parameters in both channels and the configuration with |e CT | > |e |. In the first case one should set 

e CT = e a = 0.016 ±0.003, 

to get the experimental value of the rf — > 3tt° branching fraction. In the second case we have e CT = 0.020 ± 0.004, 
e a = —0.011 ± 0.002. In fig. [5] Z-distributions of theoretical predictions of squared matrix element for presented above 
parameter values are shown. It can be seen, that these curves are different mostly at large Z, i.e. at the end of dalitz 
region, where experimental errors are large. Moreover, as it was mentioned above, in our parametrization the shape 
of Z-distribution depends only on the ratio of the model parameters, for which the uncertainties are compensated 
with good accuracy. 



Let us now proceed to tf 
written in the form 



IV. 7/ 7r°7r + 7r- 

decay. According to results of sec. II the matrix element of this decay can be 



A [77' -> 7T° {pi)ll + (j> 2 ) 7T (p 3 )] = e a Aa(s23)+e a {Aa(si2)+Aa(si3)} + e p Ap(si2,S 13 ), 

where A a , a ,p are the amplitudes of virtual cr-, a®- and p-mesons exchange and £^,a,p are isospin violation parameters 
ion these channels. Numerical values for e CTiQ were obtained in the previous section from 7/ — > 3ir° decay analysis, and 
e p can be determined from the fit of experimental data. 



1 subroutine PROB from [TJ 
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TABLE II: Coefficients Bij for rj' -> 7r 7r + 7r- decay 
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In the case of rj 7r°7r + 7r only the branching fraction is known experimentally (18j ]: 

5i? ea;p (?7' -> 7r°7r+^-) = (37±| 1 ± 4) x 10- 4 . 

It is convenient to write model parameter dependence of the theoretical prediction of this branching fraction in the 
form similar to expression ([7]): 

BR th (rf -»• 7r°7T + 7r-) = ^ <,<.,//,,. 

Zjj— <r,a,p 

Numerical values of the coefficients Bij are presented in table |TTJ Since the coefficient B pp exceeds significantly 
the others, the parameter e p should be small in comparison with e CTjQ . The possible reason is that the mixing of 
different isospin states in vector channel could occur only through three-gluon annihilation. If one neglects p-meson 
contribution, the branching fraction of rf — > 7r°7r + 7r _ decay is 

BRth (r/ -> 7r°7r + 7r-) = (22 ± 3) x 10" 4 

for e a = e e = 0.016 ± 0.003 and 

BR th {rf 7r vr + 7r") = (12 ± 2) x 10~ 4 

for e CT = 0.022 ± 0.04, e a = —0.011 ± 0.002. In the first case the theoretical prediction for the branching fraction 
agrees with the experimental value, while in the second case a large discrepancy is observed. This discrepancy 
can be removed fixing the isospin violation parameter in p-meson channel. According to our analysis, experimental 
value of the branching fraction is reproduced by setting e p = (1.6 ± 0.8) x 10 -3 in the e a = e a = 0.016 ± 0.1 and 
€ p = (2.2 ±0.2) x 10~ 3 for e a — 0.022, e a = —0.011 parametrization. Just as it was expected, in both cases e p is about 
an order of magnitude smaller than the parameter e a . One can determine what of the variants is realized in nature 
only from the analysis of the dalitz-plot shape. Unfortunately, there is no experimental information about it yet. 
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FIG. 6: V-distributions of squared matrix element of the decay rj —>■ 7r°7r + 7r~. Solid and dashed line correspond to parameter 
values e a = e a — 0.016 and e CT = 0.022, e a — —0.011 respectively. On the left figure p-meson contribution is neglected. On 
the right figure the diagram with virtual p-meson exchange is taken into account, (the isospin violation parameter is equal to 
€ p — 1.6 x 10 -3 for the solid curve and e p = 2.2 x 10 -3 for the dashed one 



The squared matrix element of this decay is usually written in the form of expansion 

\A\ 2 - 1 + aY + bY 2 +cX + dX 2 + (11) 
where dalitz variables X, Y are defined according to 

X = V3 T+ n T - , F = 3 r °~ W "° , Q = M v > - 2m„+ - m n o. 

In the above expression To and T± are the energies of neutral and charged 7r-mesons in the initial meson rest frame. 
In fig. [6] we show F-distribution of the squared matrix element for different parametrizations. Sold and dashed curves 
on this figure correspond to parameter sets e CT = e a = 0.016 and e a — 0.022, e a = —0.0011 respectively. In figure [6^ 
the contribution of p-meson was neglected, while in fig. [6Jd it was taken into account. It can be easily seen, that, in 
contrast to rf — > 37r°-decay, different parametrizations of the matrix element give significantly different distributions 
over the dalitz plot. In addition, despite of the smallness of isospin violation parameter in p-meson channel, the 
presence of this subprocess changes the distribution dramatically. For this reason experimental information on this 
distribution will be very useful. 



V. CONCLUSION 



The article is devoted to the investigation of dalitz distributions of isospin violating decays rj' — > 3tt°, 7r°7r + 7r~. 

This decays are usually studied in the framework of chiral perturbation theory. It is well known, however, that for 
high energy reactions (for example, rj' — > 37r) this theory is not valid. As a result, chiral perturbation theory does not 
describe well branching fractions and dalitz distributions for these decays [l9j . 

We think that there is another, more physically motivated, approach for dalitz-plot description. In the framework of 
so called isobar model the matrix element of the considered process is saturated by contributions of virtual resonances 
with suitable masses and quantum numbers. In our recent paper [l2[ we use isobar model to study 77' — > 777r°7r° decay. 
It is shown in this work, that dalitz plot of this decay can be reproduced by taking into account two scalar resonances: 
ao-meson in 7ny-channel and er-meson in 7T7r-channel. The integrated width of rj' — ¥ rfmr decay is determined mainly 
by ao-meson contribution, and for describing the shape of the dalitz-plot it is necessary also to take into account 
cr-meson contribution. 

In our present paper we use the same approach to study 7/ — > 37r decays. It is well known, that in these decays 
G-parity is violated, so its thorough investigation could give additional information on the nature of this effect. 
Actually, in our model there are several scales of isospin violation: either in the 77' —¥ an decay vertex or in final 
7r-meson interaction (i.e. in clq — > irir reaction). The role of each of these mechanism can be determined from the 
analysis of 77' — > Sir decay dalitz-plot. Presented in our paper results show, that coupling constant of ?/ — > an 
vertex is suppressed in comparison with isospin allowed rj' — > ar\ coupling constant by a small factor of the order 



8 



(m<2 — m u )/m s w 0.02. This value agrees well with previous theoretical estimates. As for isospin violation in virtual 
ao-meson decay, the suppression in this case is much stronger. This result is rather surprising and can be caused by 
large experimental uncertainty in the shape of the dalitz-plot. Therefore we studied also another variant, when isospin 
violation parameters are equal to each other and agree well with theoretical estimates. It is interesting to note, that 
dalitz-distributions in these cases are rather close to each other. We think that this question deserves more thorough 
experimental investigation. 

In the last part of our paper we consider the charged mode of the f] 1 — >• 37r-decay, i.e. the reaction rf — > ir°ir + ir~ . In 
this case, in contrast to neutral decay mode, p-meson exchange can also give contribution into the matrix element. The 
isospin violation parameter in this channel turns out to be about an order of magnitude smaller (~ 10~ 3 ). The possible 
reason is that the mixing of different isospin states in vector channel can occur only via three-gluon annihilation. 
Extremely interesting is the dalitz-plot shape in this decay. It turns out, that, contrary to 77' — > 3ir° reaction, 
variation of model parameters changes the form of this distribution significantly. The experimental information in 
this distribution could be very important. 

Acknowledgments 

This work was done under financial support of Russian Foundation for Basic Research (grants #10-02-00191-a, 
09-02-00132-a, 07-02-00417-a). The work of one of the authors (L.A.V.) was also supported by President grant 
(#MK-110.2008.2), grant of Russian Science Support Foundation and noncommercial foundation Dynasty. 



[1] S. Narison, Nucl. Phys. Proc. Suppl. 186, 306 (2009), 0811.0563. 
[2] N. A. Tornqvist, Z. Phys. C68, 647 (1995), hep-ph/9504372. 
[3] D. P. Majumdar, Phys.Rev.Lett. 21, 502 (1968). 
[4] J. A. Cronin, Phys. Rev. 161, 1483 (1967). 
[5] J. Schwinger, Phys. Rev. 167, 1432 (1968). 

[6] P. Di Vecchia, F. Nicodemi, R. Pettorino, and G. Veneziano, Nucl. Phys. B181, 318 (1981). 

[7] N. Beisert and B. Borasoy, Nucl. Phys. A705. 

[8] J. Bijnens, Acta Phys. Slov. 56, 305 (2006), hep-ph/0511076. 

[9] N. G. Deshpande and T. N. Truong, Phys. Rev. Lett. 41, 1579. 
[10] A. H. Fariborz and J. Schechter, Phys. Rev. D60, 034002 (1999), hep-ph/9902238. 
[11] B. Borasoy and R. Nissler, Eur. Phys. J. A26, 383 (2005), hep-ph/0510384. 
[12] S. V. Donskov, A. K. Likhoded, A. V. Luchinsky, and V. D. Samoylenko (2009), 0902.3320. 
[13] D. G. Sutherland, Phys. Lett. 23, 384 (1966). 

[14] D. J. Gross, S. B. Treiman, and F. Wilczek, Phys. Rev. D19, 2188 (1979). 

[15] B. Borasoy, U.-G. Meissner, and R. Nissler, Phys. Lett. B643, 41 (2006), hep-ph/0609010. 

[16] A. M. Blik et al., Phys. Atom. Nucl. 71, 2124 (2008). 

[17] CERN Program Library, 1998. 

[18] P. Naik et al. (CLEO), Phys. Rev. Lett. 102, 061801 (2009), arxiv:0809.2587. 
[19] J. Bijnens and K. Ghorbani, JHEP 11, 030 (2007), 0709.0230. 

[20] I. Caprini, G. Colangelo, and H. Leutwyler, Phys. Rev. Lett. 96, 132001 (2006), hep-ph/0512364. 
[21] F. J. Yndurain, R. Garcia-Martin, and J. R. Pelaez, Phys. Rev. D76, 074034 (2007), hep-ph/0701025. 
[22] I. Caprini, Phys. Rev. D77, 114019 (2008), 0804.3504. 
[23] J. R. Batley et al. (NA48/2), Eur. Phys. J. C54, 411 (2008). 

[24] T. Feldmann, P. Kroll, and B. Stech, Phys. Lett. B449, 339 (1999), hep-ph/9812269. 

Amplitudes 

1. Aa 

During last years chiral pertubation theory and Roy equations [20] make possible to desribe accurately the amplitude 
of 7T7r-meson rescattering at low energies. Unitarity relations put rather strong constrains on this amplitude. First of 
all, for s < (2m 7T ) 2 it should be real, while for s > (2vn^)' 1 up to K K prodution theshold its imaginary part should 
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satisfy the relation 



Im 



A(s) 



1 - 



4m£ 



Moreover, according to chiral perturbation theory, this amplitude should be equal to zero at s — sa — /2 (so called 
Adler effect). 

The amplitude of 7T7r-rescattering, that satisfy listed above conditions can be expressed through the series in a 
variable 



w(s) 



K 



that transforms s-variable plane with cuts s < and s > Am} K into a unit disc |wj2el. The itroduction of ap new 
variable improves the convergence of the series in the considered variable domain. The amplitude of 7T7r-scattering is 
written in the form 1211: 



Kip(s) = k ■ 



sa 



2s A 



B Q + Biw(s) + . . 



-i\ 1- 



(12) 



where unkonwn constant n can be determined from the fit of experimental data. 

The analysis of NA48/2 results in K e 4 decay [22], |23j tells us that in the presented above series one can leave only 
first two terms with the coefficients 



B = 7.4, B l = -15.1. (13) 

With these values of parameters the amplitude (fT2|) has the pole at y/s — (459 + 209i) MeV, that can be interpreted 
as a pole of a-meson with mass and width equal to 

M a = 459 MeV, r CT = 518MeV. 

It is clear, that there are other resonances (for example /o-meson) in this channel. Since the parametrization f| 12[) 
was obtained from the fit of experimental data, it includes all these contributions automatically. 

The numerical value of the k constant depends on the quantum numbers of the initial particle and can be determined 
from expeimental data. The analisys of the rf — > ri7r ir decay, presented in [liq . shows, that in rj' — > arj — > -k-kj] 
reaction reasonable agreement with experimental dalitz-plot is obseved for 

K = -4.0. (14) 

In our case there is G-parity violation in rf — > o"7r, that can be taken into account by small suppression factor. 
Recalling all mentioned above we use the followinng form of the amplitude of cr-meson subprocesses: 

e a A a = £ ff K[io(si2)+io(si 3 )+io(s23)] , (15) 

Here e CT is a small suppression factor of the order {md — m u )/m s , and = (p^i + Pirj) 2 are invariant masses of 
pairs of final pariticles. The following relation holds for these variables: 

si2 + sis + s 23 = Af2 + 3ml . (16) 



2. An 



It is well known, that the mass of ao meson is close to KK production threshold. As a result, that propagator 
of this meson has a characteristic form different from usual Breit-Wigner parametrization. One has to introduce the 
energy- dependend self-energy part, caused by nr], irrj' and KK loops, into the propagator. It is clear, however, that 
these corrections are significant only near ao-meson pole. For the allowed in 77' — > 37r-decay mass region they are not 
very important. More important is the form of aQ-r]'?: and aoirir interaction vertices. 

There are several ways to parametrise these verticies. One can use simple point-like interaction with coupling 
constants determined from experiment. Our analisys of r) 1 — > ?77r 7r -decay shows, however, that these constants do 
not satisfy the relation of 5J7(3)-symmetry [2J|. The other parametrization is connected with chiral perturbation 
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theory, that requres the interaction amplitude to be equal to zero in — ► limit. In the work [12j it was shown, 
that if one uses expressions 7^77 (PttPjj ) and j„ v i (p w p v ' ) for a irr] and aoirrj' interaction verticies, then the fit of the 
experimental dalitz-plot gives comparable values of 7-constants, as it is required by SU (3)-symmetry: 

1m, ~ 77T7)' ~ 7GeV _1 , 7^7^' = 35 GeV~ 2 . (17) 

In our present work we use the same parametrization. The only difference is that there is isospin violation in ao — > tx~k 
decay. This effect can be parametrized introducing the small suppression factor of the order (rrid — m u )/m s . So, we 
use the following expression for the a -meson exchange amplitude: 

e a A a = e 7irt77jrr)' n , . ' , r + premitations, (18) 

S12 - m a + im a T a (S12) 

Here the width of ao-meson is determined by corresponding loop integral and e a is a small parameter that will be 
determined from the fit of experimental data. 



3. A p 

In the case of charged decay mode there are diagrams with virtual p exchange (neutral yO°-meson is forbid- 
den by charged parity conservation). The vertex of p^-meson decay into it^tt pair can be written in the form 
9p (p-k+ ~ Ptt )^ where e M is the polarization vector of p-meson and g p is effective decay constant. Numerical value 
of this constant can be determined from p-meson decay width: 



T^TT ) 



24tt 



1 - 



Ami 



3/2 



Using tabular values of p-meson mass and width it is easy to obtain g p 4.2 for this constant. The vertex of 
rj' — > p^tt^ can be written in the similar form, but the corresponding coupling constant should be suppresed by a 



small isospin violation factor. As a result, the amplitude of rf 



decay with virtual p-meson has the form 



e P A p (s 13 ,s 13 ) = 2g 2 p 



S12 - S23 



812 ~ 8X3 



Sl3 



I + im p T p s 23 - m 2 + im p T p 



